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With a view of applications to the simulations of supernova explosion and proto neutron star 
cooling, we derive the Boltzmann equations for the neutrino transport with the flavor mixing based 
on the real time formalism of the nonequilibrium field theory and the gradient expansion of the Green 
function. The relativistic kinematics is properly taken into account. The advection terms are derived 
in the mean field approximation for the neutrino self-energy whiles the collision terms are obtained in 
the Born approximation. The resulting equations take the familiar form of the Boltzmann equation 
with corrections due to the mixing both in the advection part and in the collision part. These 
corrections are essentially the same as those derived by Sirera et al. for the advection terms and 
those by Raffelt et al. for the collision terms, respectively, though the formalism employed here is 
different from theirs. The derived equations will be easily implemented in numerical codes employed 
in the simulations of supernova explosions and proto neutron star cooling. 
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I. INTRODUCTION 



The neutrino transport plays an important role in some astrophysical phenomena such as supernova explosions and 
the following proto neutron star cooling (e.g. and references therein). In their studies, the Boltzmann equation 
ly-j i or its approximate versions are commonly employed to describe the temporal variations of neutrino distributions in 
{vq ■ the phase space. These equations are usually derived from the following assumptions [[^J: (1) the neutrinos are 
propagating along the geodesies for a massless particle, p^p^ — 0, and the volume in the phase space occupied by 
these neutrinos is not varied along their world line if there is no reaction. (2) the variation of the neutrino population 
due to reactions is described by the so-called collision terms obtained with the stohszahl ansatz. With the masses 
non-diagonal in the neutrino flavor space, the neutrino oscillation occurs among different flavors of neutrinos(e.g. ||] 
(*S and references therein). It is thus interesting from an academic point of view how this oscillation phenomenon is 
described by the generalized Boltzmann equations [prH] • It is also important from a practical point of view for those 



i 

2 ' who are interested in the possible significant consequences the oscillation might give in astrophysical events [[L5|- 

In collapse-driven supernova explosions, for example, this is particularly the case if the resonance of oscillation occurs 
near a neutrino sphere where neutrinos are interacting with other particles and thus the oscillation should be treated 
simultaneously with these reactions and possibly with the evolution of the matter distribution as well. The purpose 
• i-H , of this paper is to provide the formulation which can be easily implemented in those numerical simulations. 

' In considering the transport equation with the oscillation, we have to rely on a more formal derivation of the 
Boltzmann equation. This might be done in a couple of ways. Sirera and Perez |l4j], for instance, based their 
derivation on the relativistic Wigner function approach in the mean field approximation. Although they took the 
relativistic kinematics properly into account, they did not obtain the collision terms, since it is difficult to go beyond 
the mean field approximation in their formalism. Raffelt et al. |Io| , pT| ], on the other hand, obtained their transport 
equation via the density matrix approach. Although they derived the collision terms, they did not consider the 
spatially inhomogeneous system. In this paper, we derive the relativistic Boltzmann equation including corrections 
due to the oscillation both in the advection terms and the collision terms by employing the real time formalism of the 
nonequilibrium field theory [fl9|-pT|. In this approach, the dispersion relation and the collision terms are derived on 
the same basis, that is, a particular approximation for the self-energy of neutrinos, which is conveniently represented 
with Feynman diagrams. 

The paper is organized as follows. We first derive a generic form of the transport equation without specifying 
particular equations of motion of fields. Then, the formulation is applied to the neutrino flavor oscillations. In so 
doing, we ignore small corrections of the order of mj; / E v except for the terms responsible for the flavor conversion, 
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as is usually the case. Here m v and E v are typical mass and energy of neutrinos in the observer's inertial frame. 
In this limit, as shown later, the left handed neutrinos are decoupled from the right handed ones and the difference 
between Majorana mass and Dirac mass never shows up in the flavor mixing. The general relativistic corrections are 
obtained up to the leading order of A„/i?, where A„ is a typical wave length of neutrino and R is a scale height of the 
background matter distribution. 



II. FORMULATION 



A. general derivation of transport equations 



In this section we derive general transport equations for multi-component fields based on the real time formalism of 
nonequilibrium field theory by Keldysh ]l9|-pl ] . In this formalism, we introduce path-ordered products of operators 
on the closed time-path, which extends from t — — oo to t = +00 then back to t = —00. In this product, the operator 
with a time argument which comes later on the above time-path is put to the left of other operators whose time 
arguments come earlier. Accordingly the path-ordered Green function is defined as 



iG piJ {hM) = (T P ^i(ii) V#2)) 



(2.1) 



Here T p stands for the path-ordered product of the following operators. The subscript ij of the Green function denotes 
the components of the field. The bracket (• • •) represents that arguments are averaged over the ensemble specified 
by a density operator p as Tr{- ■ ■ p}, where Tr is a trace operator. We define a generating functional of the Green 
function as 



Z(J,,P) = Tr < T p 



exp(z^ J d A x (^jj (x)if)i(x) + Ji(x)tpj(x) 



exp 



[iw(j,jiy 



(2.2) 



The Green function is obtained by the functional derivative, i G P ij{x,y) = —fji 



ating functional for the connected Green function is denoted as W(J, Jt). Going to the interaction representation, 
we obtain 



Z(J, J^)=Tr\ T p 



exp 



S Jt ' iSJ 



Trl T n 



exp(i^ J d 4 x ( j/(s)^ 7i (a;) + J l {x)ip] l {3 



(2-3) 



where the Lagrangian density for interactions is denoted as Ci n t and the subscript / indicates that the variables are 
given in the interaction representation. The last factor of the right hand side of Eq. (|1^) is the generating functional 
for the no interaction case, Zq(J, J'), and is given as 



Z (J, Jt) = Z vac (J,j1) Try. exp ft j ^ ( 
with the generating functional for vacuum, 

Z vac (J, J f ) = exp I - i £ f f d A xd\j J 3 {y) G%(x,y) j}{x) 



(2.4) 



(2.5) 



Here Gp^ (x,y) is the path-ordered Green function for vacuum. The normal order product is represented by '. ■ ■ ■ '. in 
Eq. (fE4|). All the information of the ensemble is included in the last term of Eq. (2.4), N p (J, jt) = exp[z W p (J, Jt) ]. 
Its connected part (J, Jt) is in general expanded to cumulants as 

W?{J, Jt) = £ -T-i / • • • /dW • • <*VdV ■ • d A x m J{y x ) ■ ■ ■ J(y n ) j\ Xl ) ■ ■ ■ j\x m ) 

Tfi. ft. J J - 

m,n=l r 

x W? mn ( Xl ,---,x m I yi, •• -,y n ) . (2.6) 
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In the following we assume that the expansion is terminated at the quadratic order. This is true, for example, for the 
thermal equilibrium and the more general condition for this to be true can be found in the paper by Danielewicz [^2| . 
With this assumption, we can expand as usual the Green functions by the propagator which have corrections origi- 
nating from a particular ensemble. 

The Dyson equations are obtained by the Legendre transformations: 

r(^,Vl) = w(j,j^) ~ Jt-Vc - J-iPt (2.7) 

with ijj c {x) = s S j\( x ) W and tpl(x) — jjj-s W. We use the abbreviation J+ • i/j c = J2 i J^d^jj (x)tp i(x). Then the 

following relations hold: ST/Sip c (x) = —J^(x), 8Y/8ip\{x) = —J(x). The Dyson equations take the integral form on 
the closed time-path as 



d z Gp(x,z)T p (z,y) = Id z T p (x, z) G p (z,y) = S p (x - y) , (2.8) 
Jp 

where the connected Green function and the vertex function are defined as iG p (x,y) — ^t( x ) i&jty) W{J,fl) and 

T p (x,y) = <5 2 r / 5if) c (y)8ipl(x), respectively. The <5-function is extended on the closed time-path as follows: S p (x — 
y) = S(x — y) for t x and t y on the positive branch of the time-path extending from t = — oo to t = +00 and 
5 p (x — y) = —6(x — y) for t x and t y on the negative branch of the time-path that runs from t = +00 to t = —00. 
Introducing the matrix representations for the Green function and the vertex function as 



we can recast the Dyson equation in a single time representation: 

r d 4 z G c (x,z)a 3 f(z,y) = ( ' d 4 z f(x, z) a 3 G c {z,y) = a 3 5{x-y) . (2.10) 



In the above equations, the time integration runs from t — —00 to t — +00, and <r 3 = ( ^ ^ J is the Pauli matrix. 

The subscripts F and F indicate that the time arguments are both on the positive branch and on the negative branch, 
respectively, while the subscript + means that the first argument is located on the positive branch and the second on 
the negative branch, and the subscript — represents the other way around. It is clear that Gf is an ordinary Green 
function defined from the chronologically ordered product while Gp is obtained from the anti-chronological ordering. 
From these quantities, we further define the retarded, advanced and correlation functions as 

G r = G F - G+ , (2.11) 
G a = G F - GL , (2.12) 
G c = G F + Gp ■ (2.13) 

The counter parts for the vertex functions are defined in an analogous way. Using the identity IV + = T+ + T_, 
we can express T's in general as 

r± = i(B ± A) , (2.14) 

r> = D + iB , (2.15) 

I> = -D + iB , (2.16) 

r r = D + iA , (2.17) 

r a = D - iA , (2.18) 

where A, B and D are three Hermitian matrices. Solving the Dyson equations using these quantities, we obtain the 
general form of the Green functions as 

G r = [D + iA}' 1 , (2.19) 

G a = [D - iAY 1 , (2.20) 

G± = -i[D + iA^iB ± A][D - iA}- 1 , (2.21) 

G F = [D + iA^lD - iB][D - iA}- 1 , (2.22) 

Gp = - [D + iA^lD + iB}[D - iA}- 1 . (2.23) 
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The dispersion relation is obtained from D and A, while the distribution function is found from B as shown shortly. 
D, A and B can be represented in turn by the self-energy E p which is defined from the two point vertex function as 

Tp = Tpo — E p , (2.24) 

where the free vertex function is r p o(a; — y) = S(d x ) S p (x — y). Here the derivative operator is taken from the free 
Lagrangian, Cq = ip'S(d) ip. Defining again the matrix components of the self-energy in the single time representation, 
we obtain 

D = S(d x )5{x-y) - i(S F - E f ) (2.25) 
A=\i (E_ - E+) (2.26) 
B = \i (S_ + E+) . (2.27) 
The self-energy, on the other hand, is given by the relation 



£ p (af,2/) = y-i ( T p j(x) f(y) ) - S p (x - y) / 



1P7 



Sip(y) Sip^{x) 



where the currents are defined as j(x) — C-int and j^(x) = 5Ci n t/5ip{x), and the subscript 1PI means the one 

particle irreducible part. 

Now we introduce the distribution function, n. First we define another Hermitian matrix N from B as 

G c = — T^ 1 2iB T^ 1 = IV 1 iV - (2.29) 

Then it satisfies the following equation 

ND ~ DN ~ i(N A + AN) = -2iB . (2.30) 

The matrix distribution function is finally defined as 

N = 1 T 2n , (2.31) 

where the upper and lower signs are taken for Fermion and Boson, respectively. It is easily shown that this distribution 
function becomes Fermi- or Bose-distribution function in the thermal equilibrium case. In that case, n can be 
simultaneously diagonalized with D and gives the distribution functions of quasi-particles. In general, however, n 
has non-diagonal components even in the representation which diagonalizes D. These non-diagonal components are 
responsible for the flavor mixing as discussed below. Eq. ( 2.30 ) gives the equation satisfied by n: 

nD — Dn - i(nA + An) = ± i(B — A) = . (2.32) 



Using Eq. (2.26), we can rewrite the above equation, 

nD - Dn = T «)(tS+) + (tS+)(1 T n)\ - \[n^- + E_n] . (2.33) 

It is already clear that the right hand side of the above equation describes collisional processes among the quasi- 
particles. In fact, (±zE+) and (— iE_) can be interpreted as the emission and absorption rates of quasi-particle. 

The transport equation as we know it is obtained by performing the so-called gradient expansion for the above 
equation. The Wigner representation of a quantity F(x, y) is obtained by making Fourier transformation with respect 
to the relative coordinate as 

F(k,X) = jd\x-y) e ik ^F(x,y) , (2.34) 

with the center of mass coordinate X = (x + y)/2. The gradient expansion is performed by taking the Wigner 
representation of both sides of Eq. ( 2.33| ) keeping only the leading order of the derivative with respect to X. Thus, 
we obtain the transport equation as 
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dD{k,X) dn(k,X) dn(k,X) dD(k,X) 



dk„ 



dX» 



dX^ 



dk„ 



dD(k,X) dn(k,X) dn(k,X) dD(k,X) 



dk„ 



dk n 



- i [D(k,X)n(k,X) 



= -{[1 T n(k,X)][ T iV + (k,X)] 



dX» 
n(k,X)D(k,X)} 

[TiX+(k,X)][l T n(k,X)]} 



dX^ 



{n(k,X) [i£_(/c,X)] + [i£_(A;,X)] n(k,X)} 



(2.35) 



It is evident that the first row of the above equation represents ordinary advection terms while the right hand side 
stands for the collision terms. The second row, on the other hand, does not appear in the ordinary transport equation 
and we see below that this term causes the mixing among neutrino flavors. What remains now to do is to give the 
self-energy which determines not only the collision terms but also the dispersion relation, that is, D. We do this for 
the neutrino mixing in the next section. 



B. neutrino transport equation with flavor mixings 

We apply the general formulation obtained so far to the neutrino transport. The following Lagrangian density is 
considered: 

£ = j i^il^dtiipL - ^iplM„ip L - \ip L Mi%l} L L + dnt for Majorana v (2 36) 

I | ^lI^O^l + § ^n^d^R - ^ R M D tp L - ^ l MUr + C-int for Dirac v 

where the Majorana and Dirac masses are M u and M Dl respectively. The subscripts L and R stand for the spinor 
with left and right handed chirality, respectively, and ipl = Cipl with C the charge conjugation and the superscript 
T representing the transposition. The interaction Lagrangian density is denoted as £i n t- In the above equation, 
the indices for spinor components and neutrino flavors are suppressed. In the following, the flavor is denoted by the 
superscript and the spinor component by the subscript as ipf when necessary. The matrix Green functions of interest 

are ( T p vj L f ~$ L ) ) , (T p ^^)' ( T P^i^ij) and (TM^rj)- Here and in the following, i> R should be replaced by 
tpl for the Majorana neutrino. We discuss the advection part and collision part of the Boltzmann equation separately, 
since we apply different approximations to the self-energies included in them. 



1. advection part 

Following the common practice, we take the mean field approximation for the neutrino self-energy in the advection 
part, whi ch is conveniently represented by a Feynman diagram shown in Fig. |l| and comes from the second term of 



Eq. ( 2.28 ). Only scattering processes contribute to this ensemble average. In the supernova core, the scatterings on 
free nucleons, nuclei and electrons are important. The former two of them occur only via neutral currents and as a 
result, the self-energies corresponding to them are proportional to the unit matrix in the flavor space: 

£S?t =E"% t^ 7 " ( : - ^) ^] K - Kl 5 ) M - (2-37) 



V2 

sr^ab eab ^ F 



^V(l-y)£^P0*P° . (2-38) 



N 

G 



^f m = -8 ab ^^{l+^)Y,KPN{X)5 llQ , (2.39) 
X a An = Zf RL = , (2.40) 



where Eq. (2.39) is true only for the Majorana neutrino and ^p m = for the Dirac neutrino. In the above equations, 
the subscript N runs over neutron and proton, and p N stands for the nucleon number density. The similar equations 
are obtained for the scattering on nuclei. Hence, in the following, the nucleon scattering is considered. 

On the other hand, the scattering on electrons gives a non-trivial structure to the self-energy in the flavor space 
since the process occurs not only through the neutral current but also through the charged current, and the latter is 
relevant only for the electron-type neutrinos in the matter in which electrons are abundant but other charged leptons 
are not. In that case, the interaction Lagrangian density becomes 
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q% = - ^§ [ r L - r (i - 7 5 ) c e ] 7„ (r ~ r 7 5 ) ^ ] 

+ E - % ^ C 1 - V 5 ) ^ ] 7m W - 9 A 7 5 ) ^ ] ■ (2.41) 

a * 

In the above equation, g v = — 1/2 + 2 sin 2 6 W and g -4 1/2 denote the vector and axial vector coupling constants 

of the neutral current, while the charged current is also taken into account in g v = g v + 1 and g A = g A + 1. The 
Weinberg angle is referred to as 6 W here. We obtain the self-energy of neutrino in the mean field approximation as 

S« b M = W- ^7 M (1 - 7 5 ) ~9 v Pe(X)5,o 

+ S a ^S b ^ Q= r (1 - 7 5 ) g v Pe (X) S^o , (2.42) 

for the unpolarized electrons p^ , |24] |. Here the electron number density is denoted as p e . As for the other components 
of the self-energy, Ej? Ifl = T, Frl = common to both types of neutrinos, and T, Frr = —T, Fll with (1 — 7 s ) replaced 
with (1 + 7 5 ) for the Majorana neutrino and ^ Frr = for the Dirac neutrino. If the electrons are polarized in the 
magnetic field, the neutrino self-energy is modified to f25| , p6"f 

E «k = sau^^F^ (! _ 7 5) {~g-p e (X)S, a - ~ 9 A p° e (X) S, z } 

+ S a —S b — ^7 M (! " 7 5 ) {g v Pe (X)S lM - g A p° e (X)S, z } , (2.43) 

where the magnetic field is parallel to the Z-axis. The electron number density in the lowest Landau level is represented 
as It is again true that the other components of the self-energy are zero except for Sf B „ = -Ef M with (1 — 7 s ) — > 
(1 + 7 5 ) for the Majorana neutrino. It is easily understood that neutrino-neutrino scatterings can be treated just in 
the same way. 

Now that we obtain the specific form of the neutrino self-energy, we can apply it to the left hand side of Eq. (2.35). 
Sup pressi ng the flavo r and spinor indices and writing only the chirality components in matrix form, we obtain D in 
Eq. (H|) using Eq. ( ^25| ) as 

( D LL D LR \ _ ( k,r $7° - $ ul z -Mt \ 

u ~ \d rl d rr J ~ \ —M + $' 7 ° + v B y- j ■ [ZM > 

Here the potentials are defined as 

$= S av 'S hv ' V2G F g v Pe + 5 au ^5 bu ^V2G F g v p e , (2.45) 
$ B = -6 av 'S bv °V2G F g A p° e - S av ^5 hv ^V2G F g A p° e , (2.46) 

with <£>' = $ and $' B = $ B for the Majorana neutrino, and $' = and $' B = for the Dirac neutrino. It is understood 
in the above equations that $ B — in the case of no magnetic field. The dispersion relations for quasi-particles are 
obtained from the eigen values of D. 

We first make an order estimate of each term in the advection part. Defining R as a typical length scale of the 
matter distribution, the density scale height, for example, E v as a typical energy of neutrino, and Am 2 as a square 
mass difference, we find 

(2.47) 

n 



dD dn 


n 




'dk'dX^ 


" R ' 




dD dn 


$ n 


Ami 


dx'dk 


R E u 


E 2 

V 


D-n r 


Am 2 


R n 
J ' R 



R ■ (2-48) 

(2.49) 

Here A is a wave length corresponding to Am^jE v : A ~ 0.1cm[Am 2 /leV 2 ]~ 1 [£ ; ^/lMeV]. For the typical mas s 
difference and energy of neutrino, Aml/E 2 ~ 10 -12 . Hence the second term in the left hand side of Eq. (2.35), 
which represents the potential force exerted on neutrino by surrounding matter, is much smaller than the first term, 



G 



which corresponds to the ordinary advection term in the Boltzmann equation. We ignore the former in the following 
discussion. 

Next we show that the n LL can be decoupled from the other components assuming that Am^/^ is neglected. We 
perform two matrix manipulations for D ■ n — n ■ D : (1) multiply the first row with /c M 7 M + $'7° + <&^7 Z from 
the left and add to it the second row multiplied with M' from the left. (2) then multiply the first column with 
&/i7 M + $'7° + $' B 7 Z from the right and add to it the second column multiplied with M from the right. Taking into 
account that (Am u / E u ) n LL and <&/E u ~ Q> B /E U ~ Am 2 /E 2 , we obtain the equation for n LL as 

- M 2 - hrf ($7° + $b7 2 ) + + $ b7 z ) k M -f] n LL k^ 



- k^n LL [W — M 2 + ($' 7 ° + $' fl7 z ) - ($7° + $ b7 z ) k^\ 

The same manipulations are done for dD/dk^ ■ dn/dX^ + dn/dX^ ■ dD/dk^ to obtain 



k^ 7" 



dn LL 
dX" 



fc<r7 <T + kaj a 



dX» 



(2.50) 



(2.51) 



In the next subsection, it is shown that n LL can be also separated from the other components in the collision terms 
by applying the same procedures. 

To the leading order of Am u /E v , Ti is a scalar with respect to the spinor index. The familiar form of advection 
terms in th e Bol tzma nn equation is obtained by taking the trace with respect to the spinor indices after multiplying 
Eqs. (|2.50D and ( |2.5l| ) with 7°(1 - 7 5 )/2 from the left : 



ft" 



d ) 



dX* 



(M 2 



V 2fc° 



L.Z 

fc° B 



M 



2cb 



2fc° 



k° E 



(2.52) 



Here the indices of flavor are explicitly included. From this equation we see that the resultant equation is identical 
for the Dirac and Majorana neutrinos up to the leading order of Am v /E v and that D is effectively replaced by D e g 
in the flavor space of the left handed neutrinos : 



D, 



1 



cff 



( k^k^ 



M 



2fc°$ 



2fc z $ f 



(2.53) 



The positive and negative zeros of D e g correspond to the energies of the neutrino and the anti-neutrino, respectively 
The transport equation for the anti- neutrino is obtained with the replacements : fc M —> — k^ , 

till * t^ll hi Eq. ( |2.52[ ) . 

In the following, we consider only the transport of on-shell neutrinos. Since we are not interested in the small difference 
~ Am 2 / E v of the o n-shel l energies among different flavors except in the terms responsible for the flavor mixing, we 
take k° = \k\ in Eq. 

In order to illuminate the structure of the advection part of t he tra nsport equation obtained above, we discuss only 
the two-flavor case of electron- and muon- neutrinos. Then Eq. ( 2.52| ) multiplied with i becomes on the flavor basis 



d 1 



dX» 



k°[H, 



with 



H = — 



AM 2 



k° E 



1 

2fc° 



1 — cos 26>o — sin 20 o 

Ao . „„ A 

\ — sin 26>n cos 20n 

\ 2 2 



V2G 



FPe 




(2.54) 



(2.55) 



where AM 2 is a mass matrix with a diagonal matrix (mf + m|)/2 • 1 subtracted. The eigen values of M 2 are m\ and 



TO2 with the latter larger than the former, and their difference is defined to be Ao 
vacuum is denoted as 9q and the following relation holds : 



M 



U 







u 



COS Uq 

sin #0 



sin c/o 
cos Oq 



mf. The mixing angle in 

(2.56) 
(2.57) 



Eqs. (2.54) and (2.55) are the same equations as obtained by Sirera and Perez p4| 
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Taking the bases on which H in Eq. (2.55) is diagonalized at each point in space, 



Ut(x) H(x) U M (x) = - 



M 2 (x) 
2fc° = 



1 

2fc" 



M?(x) 
M 2 (x) 



with 



we obtain 



Ut(x) 



U M (x) 



U M (x) 



cos 9 M (x) 
sin 9 M {x) 



sin 9 M (x) 
cos 9 M (x) 



d i 



dh 



h, Ul{x)ik» 



dU M (x) 



ik» 



d9 M (x) 
dX» 



(2.58) 

(2.59) 

(2.60) 
(2.61) 



Here the mass eigen values and mixing angle in matter are denoted as Mi, Mi and 9 M , respectively, and the distribution 
function in this representation is defined as n. In order to see the oscillation among different flavors, we write down 
each component of the above equation : 



ik» 



dh 11 


- ik" 


d9 M (x) 


dX^ 


dX» 


On 12 


+ ik^ 


d9 M {x) 


8X» 


dX» 


dh 21 


- ik^ 


d9 M {x) 


dX^ 




dX» 


On 22 


+ ik^ 


d9 M (x) 


8X» 




dX» 



^12 



( n 
(n 12 



-21 

n 



-22 ' 

n 



-21 

n 



Am 






h 


2 




A M 






h 


2 





(2.62) 
(2.63) 
(2.64) 
(2.65) 



where the mass square difference in matter is defined as A M = M| — M 2 . Ignoring the collision terms for a moment 
and adding Eqs. (2.62) and (2.65), we obtain the relation 



d 



dX» 



n 11 



-22 ' 

n 



= 



(2.66) 



which expresses the number conservation of neutrinos. From Eqs. ( 2.63| ) and (2.64) the following equations are 
obtained, 



d 



dX» 
d 

dX» 



A, 
2 

A, 
2 



+ ik 



■ 9 9 M {x) , _., | . 



dX^ 



n — n 



= 



Although we can infer the oscillating nature of the solution, it is better seen by eliminating n 
only the leading terms of Am v jE v . The resultant equation roughly becomes 

,2 / \ \ 2 



^21 



dt 2 



■ - 12 

n 



-21 

n 



— ( 



-21 ' 

n 



1 



-22 ' 

n 



(2.67) 
(2.68) 
and taking 

(2.69) 
+ h 21 have 

an oscillating part with an oscillation length of A -1 = A M / 2E V and a non-oscillating part which is negligible when 
the adiabatic condition X/R <C 1 is fulfilled. Hence we can ignore the non-diagonal components n 12 and n 21 of the 
matrix distribution function if we are interested only in the variation of the neutrino population on the length scale 
much longer than A and consider only the mass difference and energy of neutrino which satisfy the above adiabatic 
condition, as is usually true for the supernova cores and proto neutron stars. The non-diagonal components n 12 and 
n 21 can be ignored in the collision terms after taking the average of the rapidly oscillating terms over the length scale 
much larger than A. 



where I is the path length and R is the typical scale length of matter distribution. It is evident that n 1 
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In the following, we set n 12 = n 21 = and consider the equations governing the diagonal components of the matrix 
distribution function for neutrinos. Following Raffelt et al. fll|] ; we represent n 11 and n 22 in terms of n v " and n^^, 
the diagonal components on the flavor basis. From the relation 



Ut n U M 



we obtain the distribution functions on the flavor basis as 







n u ' ±tan26> M (n^ - n v » ) 

" ~ 1 | tan 26 u {n v ° - n"" ) 



Inversely transforming the equation for n in this approximation, 

dn 



(collision terms) , 



(2.70) 



(2.71) 



(2.72) 



we finally obtain the equation for n in the limit of the adiabatic oscillation between two flavors as 



dn 
dXt 1 



+ jtan20 M {n v ° - n v » ) ( J ^ j - (n 1 " - n v " ) ( J 



nl 

OH** 3Xm 



(collision terms) 



(2.73) 



Here the mixing angle in matter is given as tan 29 M = sin 29 / { cos 29 — [ 2\[2 Gf ( p e — k z /k° p° e ) E v ]/Ao }. 

We briefly discuss here the general relativistic correction terms in the advection part. We take an arbitrary 
point in space-time and consider a small patch of space-time around it of the size of t which satisfies the condition 
R ^> £ S> 1/E V . Then we can take a local coordinates X R in this small region, which has a Minkowskian metric up to 
the second order of i/R. We also define an orthonormal tetrad e l R aligned to this coordinate and use it to project the 
four momentum of neutrino on it. On this coordinate in the small patch of space time, the above derivation for the 
advection terms in the Boltzmann equation is still valid, that is, we obtain Eq. ( p. 52] ) or Eq. ( p. 73 ) for the adiabatic 
two-flavor mixing with trivial replacements of X 11 with X% and fc M with k%. Thus, there is no additional mixing due 
to the general relativistic gravity under the current assumpti on th at the ti ny m ass difference of neutrinos is ignored 
except for the mixing term, that is, the second terms of Eq. (2.52) or Eq. (2.73). All we have to do now is to make 
a coordinate transformation X% — > X 11 and an associated momentum transformation fc^ — > fc M . The latter is, in fact, 
induced by the transformation of the tetrads, e l n — * e l , the latter of which is given globally. Employing the orthogonal 
transformation between two tetrads e % {X) = T*(X) e J R and the transformation matrix 



dX» dW 



dX v R dX v R 
dXf d W 



/ dX^ 
dX» 



dXP dkl 



dX» dXP y 

J V (T-i)J 
we can perform the transformation for the advection term as follows: 



T 



-i \3 



\ 



(2.74) 



dn 
dXj 



k % 



k l e 



dn 
dX» 
dn 



k l k' 



d< 



dXP 

n dn 

dk n 



dn 
dk^ 



(2.75) 



In the above equation, the inner product of two vectors is denoted as (v\ ■ V2) and the component of the connection 
1-form is de signated as uifj — (V ei ej ■ e k ) with W ei the covariant derivative in the direction of ej. The second term 
of Eq. ( g7F§ ) is a familiar correction term due to the general relativity [||-^) , which accounts for the red shift and ray 
bending of neutrino in the gravitational field. Since the mixing term and the collision terms (see b elow ) do not contain 
spatial derivatives, they are unaffected by the above transformation. However, d9 M /dX R in Eq. ( [2.73 ) is affected just 
in the same way as dn/dX R shown above. It is noted that this term actually originates from the advection term 
dn/dX R due to our pointwise choosing of the local mass eigen state basis. What remains to be done is, thus, to 
calculate the connection 1-form. 
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2. collision part 



In this section, we derive collision terms in the Born approximation for the neutrino self-energy. It is well known 
that the approximation of the self-energy for the advection terms is different from that for the collision terms [p7| . 
The Born approximation is conv eniently represented by the Feynman diagram shown in Fig. |^. Only the first term of 
the right hand side of Eq. (2.28) contributes to £± in the collision part. As done in the previous section, we evaluate 
the self-energy coming from various processes separately in the following. 

For the nucleon scattering, the self-energies are obtained as 



iT^Jx,y) = - 7 5 )^-^,y)f(l - 7 5 )] S mv {x,y) , 

iXf LL (x,y) = ^ [ 7 " (1 - 7 5 )* Gf LL (x,y)Y(l ~ 7 5 ) ] S N ^(y, x) . 



(2.76) 
(2.77) 



In the above equations, the dynamical structure function for nucleon is defined with the weak neutral current of nucleon 
Jn(x) as S^(x,y) = (J^(x)J^(y)). The weak neutral current for nucleon is given by J% = "0 N 7 M {h^ — h^J 5 ) ip N - 
The other components of the matrix Green function are zero for the Dirac neutrinos. For the Majorana neutrinos, 
E LK is obtained, for example, by replacing G LL with G LR and (l — 7 5 ) with 7" (l + 7 s ). Recalling the relations 
n LR ~ n RL ~ (Am v /E v )n LL ~ (Am 1/ /E u )n RH and £ iJ? ~ S fii ~ (Am v /E U )T, I , L — (Am„/£„)S Bii , we find that the 
LL-componcnt can be decoupled from the other components after the same matrix manipulations as done for the 
advection part and that the resulting collision terms are identical for the Dirac and Majorana neutrinos, if we take 
only the leading terms of Am v jE v . Note that the exchanged terms are added in the collision part unlike in the 
advection part. From the term n ■ for example, we obtain [fc^7 M n LL ] ■ [i'E^. LL k^j tJ/ ] as the LL-component. 

Following the procedures taken for the advection part, we multiply the collision terms with 7 (l — 7 s ) /2 from the 
left, take the trace with respect to the spinor indices and divide by 4fc°. We then obtain from n ■ for example, 
the following: 



(0i i ^7 1 r|7°^^fc P 7 P n 11 ^ 7 ^l-7 5 ) l G- LL {k-q,X)Y(l-l 5 ) k^] S NIXU {q,X) 

Ignoring again the tiny masses of neutrinos and A in deriving G± from D, A and n, we obtain 
1 - 7 5 



iG ±LL {k) 



krf 2ir S(k 2 ) 



(2.78) 



(2.79) 



In the above equation, it is explicitly indicated that the negative energy contribution to the number density of the 
neutrino corresponds to the number density of the anti-neutrino. It is noted that the number density is a function 
of k after we ignored A and imposed an on-shcll condition k 2 = 0. The upper (lower) components in the columns 
correspond to G+ (G_). Inserting this relation to Eq. (2.78) and recalling that n is a scalar with respect to the spinor 
indices, we obtain the following collision term for the neutrino distribution function: 



d 3 k' 



1 1 



n y (k) [ 1 - n u (k') 



G 



L^(k,k')S mu (q,X) 



(2tt) 3 2k'° 2 1- — v- /j 2 

where k! = k — q is the four momentum of the scattered neutrino, and the tensor L pv is given as 
L^(k,k') = Tr{k p j p Y l (l ~ l 5 )Ki a j u {l - 7 5 )} = 8 |Vfc v + k v k p - g^k p k' p - ie^kpk'A 



(2.80) 



(2.81) 



Here the metric tensor is denoted as g pv and the anti-s ymme tric tensor as £^ v f a with e 0123 = 1. From the term • n 
we obtain the collision term which is obtained from Eq. (2.80) by replacing n v {k) [ 1 — n„(fc') ] with [ 1 — n v (k') ] n„(fc). 
Just in the same way we obtain from the term [1 — n][— the following collision term: 



1 1 



d 3 k' 



(2tt) 3 2k'° 2 



[1 - n„(k) ] n u (k') % L""(fc, k') S mii (-q,X) 



(2.82) 



For the term [— i£+][l — n] we replace [1 — n v (k) } n v {k') with n u (k') [ 1 — n„(fc) ] in the above equation. Using the 
relation S Nufi (—q) = e~@ q S Npv {q) for the matter in equilibrium, which stands for the detailed balance, we finally 
obtain the collision terms as 
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' <Pk' 1 1 [ _ p j> n v {k')[l - n v {k)] + [1 - n v {k)]n v {k') 
(2tt) 3 2k'° 2 \ g 2 



[ 1 - n v (k') ] + [ 1 - n v {k') } n v {k) \ G% 

— (fc, k ) S mu {q,X) 



(2.83) 



If the matrix distribution function is diagonal, the above equation reduces to the ordinary collision term. 

If the mixing occurs adiabatically, the above collision term is further s impli fied. For t he tw o-flavor case (y e and 
Vp, for example), we insert the matrix distribution function given by Eq. (2.71) into Eq. ( 2.83|) . Then we obtain the 
collision term for the v e distribution function as 



1 1 



d 3 k' 



(2tt) 3 2k'° 2 



-0q° 



n v '{k')[l-n v °{k)] 



n v '(k)[l-n u '(k')[ 



i tan26» M (fc') [n Ue (k') - n v »{k')] ^ tan26> M (fc) [ n u " (k) - n v »(k)] 
- tan 26 u {k) [ n v " (k) - n u " (fc) ] - tan 26 M (k') [ n v ' (k') - n"" (fc') ; 



x ^ L^{k,k')S mv {q,X) 



(2.84) 



The collision term for the f M distribution function is obtained by replacing n v " with n Vl " in the above equation. It 
is noted that the correction terms due to mixing cancel each other for iso-energetic scatterings, which we commonly 
assume for the neutrino-nucleon scattering in the supernova cores and proto neutron stars P, |2q , ^9| . 

The collision terms for the neutrino-electron scattering are essentially the same as those obtained for the nucleon 
scattering. The main difference originates from the fact that the electron weak current has flavor dependence, which 
gives rise to non-trivial contractions of flavor indices between the electron structure function and the neutrino distri- 
bution function such as n" c (fc) [l — S^, where the superscripts a, 6, c represent flavors. The structure function S e 
is an electron counter part of the nucleon structure function S N and is defined as S^ v (x, y) = (^(^-^(y)) ■ Here the 
weak current for electron is given as Jj? = ip e 7 M (g v — g A 7 5 ) ip e for the v e scattering and J£ = ip e 7 M [g v — g A 7 5 ) ip e 
for the and v T scatterings, respectively. As a result, the collision term for the neutrino-electron scattering becomes 
for the electron type neutrino in the case of the adiabatic two-flavor mixing as 



d 3 k' 1 1 



(2tt) 3 2k'° 2 



-^V'(fc')[l - n«*(fc)] - ««*(*) [1 - n"°(k')}] S%r(q,X) 
_ i 



- taxi2d u (k')[n u '(k') - n u »(k')] ^ tan20 M (fc) [n Ve {k) - n^(k)' 



Seffi(q,X) + SZu"{q,X) 



L^{k,k') 



(2.85) 



The Vn counterpart is obtained by the replacement of v c 



in the flavor indices in the above equation. 



Although we assumed in the above derivation that the four momentum transfer is space-like to describe the scat- 
terings, it is obvious that the same Feynman diagram represents the annihilation and creation of neutrino pairs if the 
transferred four momentum is time-like. As stated above, the transport equations for the anti-neutrinos are obtained 
from the negative energy part of the distribution function, in that case , it is noted that the mixing angle 6 U should 
be also calculated for the negative energy. As is obvious from Eq. ( 2.55 ) the sign of the potentials is changed for the 
anti-neutrinos and the resonance conversion does not occur in this case as is well known. It is noted that neutrino- 
neutrino scatterings are treated just in the same way by substituting the neutrino structure function, which in turn 
should be evaluated with the neutrino Green functions, Eq. ( [2.7S ). 

Next we consider the neutrino emission and absorption on nuclcons. For the temperature and neutrino energy 
of current interest, the muon is not abundant and only the electron-type neutrino is involved in this process. The 
interaction Lagrangian density is 



fj.s 



rabs 
Lj int 



"72 



~i> c ^ (1 - 7 5 ) w ] [ i> P 7„ (g v - g A i b )^] + h. c. 



(2.86) 



where the coupling constants g v — 1.0 and g A 
approximation, the self-energy is given by 



1.23, and H. C. stands for the Hermite conjugate. In the Born 
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i^°(x,y) = ^ [ 7 M (1 - 7 5 )* £1(^)7^(1 - 7 5 )] , (2.87) 

iT%2(x,y) = j 5 )iG e + (x,y)~f»(l - 7 5 )] S np ^(y,x) , (2.88) 

where the Green functions for electrons are denoted as G§_ and the structure functions for the charged weak 
currents of nucleons are defined, for example, as Sg%(x,y) — (Jg n (x) J^piv)) with the charged current given by 

Jnp = typlii \9 V ~ 9 A 7 5 ) V'n- The Green functions G± are essentially the same as G± LL in Eq.( 2.79 ) with the self- 
evident substitution of n v with the electron distribution function f e . Following the same procedure as shown above 
and employing the detailed balance relation satisfied by nucleons in thermal equilibrium, 

SZi-l) = e-««° + A ^ S£(q) , (2.89) 

with the difference of the chemical potentials, Afi np = fj, n — [i p , we obtain the collision term for the emission and 
absorption of neutrinos on nucleons as 

W e l{ e ~ 0[q ° + AtlnP) feiPe) [1 " - n " {k) [1 " MPe)}}^L, u (k,p e )S^(q) , (2.90) 

in the adiabatic mixing case. Here p e and E e are the momentum and energy of electrons, respectively, and the 
transfer four momentum is q = k — p e . Since there is no other components of self-energy in flavor space than the v e v e 
component, the resulting term is identical to those with no neutrino mixing. 



III. SUMMARY 



With a view of application to the simulations of supernova explosion and proto neutron star cooling, we have 
derived a Boltzmann equation with the neutrino flavor mixing being taken into account. The derivation is based 
on the nonequilibrium field theory, and the ordinary gradient expansion has been performed. We assumed that the 
typical neutrino wave length is much shorter than the scale height of the background matter distribution, which is 
true for the supernova cores and proto neutron stars. The neutrino distribution matrix which is non-diagonal in the 
neutrino flavor space is introduced. Following the common practice, the advection part has been obtained in the mean 
field approximation, where the self-energy of neutrino is non-diagonal in the flavor space. This self-energy gives rise 
to the term in the advection part, which is responsible for the neutrino mixing and does not appear in the ordinary 
transport equation. The collision terms, on the other hand, have been calculated in the Born approximation. The 
collision terms also have corrections due to the mixing. In these derivations, the relativistic kinematics is taken into 
consideration. We have further simplified the Boltzmann equation for the adiabatic flavor mixing, which is a good 
approximation in the supernova cores and proto neutron stars. The advection terms thus derived are essentially the 
same as those derived by Sirera and Perez JT^| , although they employed the Wigner function formalism in the mean 
field approximation and did not give collision terms. The collision terms derived here, on the other hand, have the 
same structure as those found by Raffelt et al. jH],[ll| in the non-relativistic density matrix method. We have also 
shown the general relativistic correction term which accounts for the red shift and ray bending in the gravitational 
field and is commonly taken into account in the supernova and proto neutron star simulations. 

The applications of the Boltzmann equation found here remain to be done. Since the corrections due to the flavor 
mixing are rather minor, particularly in the case of the adiabatic mixing, it will be simple to implement them in the 
neutrino transport code we have now at our disposal Q. This is already underway. Since the mixing angle in matter 
is dependent on the neutrino energy and the direction of momentum with respect to the magnetic field if it exists. In 
the analyses of the neutrino flavor mixing in the supernova core, it is usually assumed that the neutrinos are flowing 
out radially |l7],[l8|]. However, they have an angular distribution near the neutrino sphere. Different positions of the 
resonant conversion due to different directions of flight of neutrinos will lead to the reduction of the neutrino flavor 
conversion. This will also be true in the absence of the magnetic field if the energy distribution of neutrinos and the 
coupling between neutrinos with different energies are taken into account. These possibilities and their implications 
to the mechanism of the supernova explosion, kick velocity of pulsars, and nucleosynthesis of heavy elements will be 
studied in the forthcoming papers. 
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FIG. 1. The Feynman diagram for the self-energy of neutrino in the mean field approximation. The dashed line stands for 
the weak interaction and the thick line represents the Green function for the particles indicated in the figure. 
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FIG. 2. The Feynman diagram for the self-energy of neutrino in the Born approximation. The dashed line stands for the weak 
interaction and the thick lines represent the Green functions for the particles indicated in the figure. The part with a shaded 
region denotes the structure functions (see the text for definitions) symbolically. 
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